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Abstract. The correlation functions near higher-order glass-transition singularities are discussed for a schematic two- 
component model within the mode-coupling theory for ideal glass-transitions. The correlators decay in leading order like 
— ln(f/T) and the leading correction introduces characteristic convex and concave patterns in the decay curves. The time scale 
T follows a Vogel-Fulcher type law close to the higher-order singularities. 



INTRODUCTION 

Mode-coupling theory for ideal glass-transitions 
(MCT) describes the transition from a liquid to a 
glass as a bifurcation in the equation for the long- 
time limit of the density-autocorrelation function 
= (p^(f)P<?>/(|P<7p) for density fluctuations of 
wave-number modulus q = \q\. If the long-time limit 

= lim(^oo 0(/(r) reaches zero, = 0, the system is 
in the liquid state. For > 0, the system is in a glass 
state [1]. The theory has been worked out in detail for 
the hard-sphere system (HSS) [1, 2] where the results of 
experiments [3] and computer simulations [4] support 
the validity of MCT. In addition to liquid-glass transi- 
tions, MCT also allows for bifurcations of higher order 
[5]. It can be shown that only bifurcations of a certain 
hierarchy can arise from MCT [6] which are equiva- 
lent to the singularities in the real roots of polynomials 
of order £ upon variation of parameters [7]. The simplest 
singularity A2 or fold generically occurs when only a 
single control parameter is changed and the liquid-glass 
transition is identified with this bifurcation. At the crit- 
ical value of the control parameter the long-time limit 
fq jumps from zero to f^'^. Close to the singularity and 
for small — fq\, the correlation functions can be 

expanded in asymptotic series yielding power laws. A 
detailed analysis including leading and next-to-leading 
order results was given for the HSS [2, 8, 9]. 

An A3-singularity or cusp can appear when two con- 
trol parameters are varied. The A3-singularity is the end- 
point of a line of A2-singularities which are identified 
as glass-glass-transition points where a first glass state, 
characterized by > 0, transforms discontinuously into 
a second glass state with /J' > . At the A3, this discon- 
tinuity vanishes and the glass-glass-transition line ends. 
Such cusp was predicted for particles interacting with a 
hard core and a short-ranged attraction [10, 11] where 



density and strength of the attraction are the control pa- 
rameters. The glass-glass transition takes place between 
a state dominated by repulsion as in the HSS and a state 
dominated by attraction, the latter arrested state was pro- 
posed to be related to a gel [11]. 

Tuning a third control-parameter, the extension of the 
glass-glass-transition line can be varied and once the 
latter vanishes, the A3-singularity merges with the liquid- 
glass-transition line and gives birth to a swallowtail or 
A4-singularity. If the range of the attraction is of the 
order of 5% of the hard-core diameter such singularity is 
predicted for a square-well system with a strength of the 
attraction of several T and a density which is slightly 
larger than the one for the glassy arrest in the HSS [12]. 

In contrast to the power-laws at the A2-singularity, the 
dynamics close to higher-order singularities is ruled by 
the logarithm of the time [5]. The wave-vector depen- 
dent asymptotic solutions along with the corrections have 
been calculated in full generality [13]. The asymptotic 
laws were demonstrated for the A3-singularity of a two- 
component model [13], and for the higher-order singular- 
ities of microscopic models with short-ranged attraction 
[14]. Logarithmic decays compatible with the predicted 
scenarios were found in recent computer simulation stud- 
ies [15, 16]. 

In the vicinity of an A3 -singularity, different liquid- 
glass-transition lines cross, and the dynamics in the liq- 
uid regime close to this crossing is influenced by three 
different singularities. Dynamical scenarios predicted by 
theory in that region [10, 12, 17] were found in experi- 
ments [18, 19, 20, 21, 22, 23] and computer-simulation 
studies [19, 24]. These findings indicate that such rich 
dynamics is indeed relevant for these colloidal systems 
and that asymptotic expansions at the MCT-singularities 
can interpret the scenario qualitatively and explain the 
data even quantitatively [17]. 

In the following, the two-component schematic model 



used already in Ref. [13] shall be reconsidered with 
respect to the A4-singularity and the evolution of the 
time scales. The advantage of the schematic model is 
that all transition points can be calculated analytically. 
Even in this simple model typical features of the full 
microscopic theory are identified. First, the asymptotic 
laws for the higher-order singularities and the schematic 
model shall be introduced briefly, further details are 
found in Ref. [13]. Second, the asymptotic approxima- 
tion is shown to work well for the correlation functions 
at the A4-singularity. Third, the time scales for the log- 
arithmic decays at A3- and A4-singularity are discussed 
for specific paths in control-parameter space. 



ASYIVIPTOTIC LAWS 

The A2-singularities are characterized by a two-step re- 
laxation around a plateau value, and the asymptotic de- 
cay laws are given by scaling functions [2]. For short 
times the decay onto the plateau is given by a critical 
decay <l>q{t) — fq{t) = hq{to/tY, while for long times one 
gets the von Schweidler law <^q{t) — fq{t) = —hq{t/t'„)^. 
The exponents are determined by a single exponent pa- 
rameter A = r(l - af/T{l - 2a) = T{1 + bf /T{1 + 
2b) < 1. For A = 1, the asymptotic solution by power 
laws becomes invalid and a higher-order singularity oc- 
curs. The dynamics at higher-order singularities A/, / > 2 
is described up to next-to-leading order by [13] 

M') = {fq+fci)+hq[(-B+B,)\n{t/x) 

+ {B2 + K^B^)ln^{t/T) (1) 
+B3\n^(t/T)+B4\n^(t/T)] . 

The wave-vector-dependent numbers f^, hq, and Kq are 
characteristic for a specific singularity, while B, Bi, and 
fq are in addition functions of the separation parameters 
El and £2 quantifying the distance of the actual state from 
the higher-order singularity [13]. The leading order in 
Eq. (1) is given by the first line with B] = fq = and the 
leading order prefactor for the logarithm is, B ^ VTsiI- 
The time scale t is fixed by matching the asymptotic 
approximation with the numerical solution for ^q{t) at 
the plateau {f^+fq). 



THE MODEL 

To mimic the (^-dependence, which is an important fea- 
ture of the asymptotic expansions above, we use a two- 
component model that was introduced for the description 
of a symmetric molten salt [25]. The model has three 
control parameters which are combined to the vector 
V = (vi, V2, V3). We will use Brownian dynamics, so the 



equations of motion for the correlators (j)q{t), q — 1,2, 
read 

T:qMq{t) + %{t)+ I mq(t-t')d,,<^q{t')dt' = G, (2a) 

mi{t) = vi^l{t)+V2^2{t), (2b) 
m2(0=V3«|>i (002(0- (2c) 

For the long-time limit of Eq (2), one gets a parameter- 
ized representation of the transition surface [13], 



fi=y- 



3-{2 + x)y 
2(l-y)2y(2-xy)' 
x^y^y^ — 2y^) 
2(l-);)2(jcV-3xy-F2) " 



(3a) 

(3b) 
(3c) 



The variables x and j with x > 4 and 1 /2 < y < 3/ (2 -Fx) 
serve as surface parameters. The exponent parameter 
A = 1 — jLl2 is determined by 



^2 



(2x2 + Ax)y'^ - 12x3^2 + (2x + A)y 



(3d) 
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FIGURE 1. Glass-transition diagram for the two-component 
model for the cut V3 = V3. Liquid-glass transition points are 
shown by the full line, the A4- singularity is shown by a star 
The vanishing separation parameters ei and £2 are shown by 
dashed and dotted straight lines, Chain lines mark the location 
of vanishing correction B2{q) = B2 +KqB^ foiq= 1,2 with 
paths indicated by crosses and pluses, respectively. 

The A3-singularities are given by A = 1 or equivalently 

jtl2 = in Eq. (3d). For x large enough, such a solution 
exists with two roots, yi{x) < y2{x), where only y2{x) is 
relevant. For small x, no such solution exists. Varying 
X, the two cusp values 3'i(x) and >'2(x) coalesce at x* 
with yi{x*) = yiix*) — y*, defining the A4-singularity, 
X* = 24.779392 ...,/= 0.24266325 .... 



The cut through the transition surface for V3 = x* is 
shown in Fig. 1 as pair of light full lines joining at the 
A4- singularity which is indicated by a star, (vj, Vj, Vj) = 
(3.132,0.1957,24.78). Attached to the A4-singularity 
we find the lines of vanishing separation parameters, 
ei (V) = (dashed) and e2(V) = (dotted), which repre- 
sent a local coordinate system. The correction ampUtudes 
at the A4-singularity are Ki = 0.3244 and K2 = -2.109; 
these yield two lines of vanishing quadratic correction 
shown by the chain lines in Fig. 1 . 
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FIGURE 2. Logarithmic decay at the A4-singularity on the 
path (x) with n = 1, 2, 3 in Fig. 1 with vi - = l/2"+l, 
v| - V2 = 0.25009/2"+ 1. The solutions of Eq. (2) are shown 
as full lines, the approximation (1) as dashed lines. Long hori- 
zontal lines exhibit the critical plateau values /*, short lines the 

corrected plateau values +fq. 



ASYMPTOTIC APPROXIMATIONS 

We first analyze the path labeled n (x) in Fig. 1. The so- 
lutions of Eq. (2) are shown as full curves in Fig. 2. The 
approximations (1) for the correlators are displayed as 
dashed lines. The time scale T was matched for (j)i (f ) and 
^2(0 independently at the corrected plateau values. The 
approximation describes the decay around the plateaus 
/* reasonably for both correlators. Since B2 + K^B^ is 
set to zero on the present path and £3 = B4 = for an 
A4-singularity [13], only the first line in Eq. (1) is rele- 
vant for the upper panel. This approximation describes 
the decay of (^i(f) from 0.4 down to 0.1 and for a win- 



dow increasing window in time with increasing n. Since 
K2 < K\, the prefactor for the quadratic correction in 
Eq. (1) for (^(f) is negative and hence the decay around 
the plateau /| is concave which is described well by the 
asymptotic approximation. However, for both q the pref- 
actor {B —Bi) slightly overestimates the absolute slope 
of the solutions. This is due to a relatively large positive 
next-to-leading-order correction Ci that renormalizes the 
prefactor to (B - B 1 - Ci ) [26] . 

Another feature of the curves in Fig. 2 is noteworthy. 
There appears a window in time outside the transient dy- 
namics between f « 1 and t ~ where the description 
by Eq. ( 1 ) is not applicable. This is caused by the close- 
by A2-singularities on the chosen path, cf. Fig. 1, where 
fq> 7 ? = 1 ) 2 which introduces an additional slow- 
ing down of the dynamics before the logarithmic decay 
is encountered. In addition, we observe that (^2(0 varies 
almost Unearly in In? between ? » 1 and t w 250. It is, 
however, easy to distinguish that effective logarithmic 
variation, that stays the same upon further approaching 
the singularity, from the characteristic behavior of the 
correlators around the plateau where the prefactor of the 
logarithmic decay vanishes with the square root of the 
distance. The almost linear relaxation in log? seen for 
^2(0 is related to the /3-peak phenomenon [27]. 
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FIGURE 3. Logarithmic decay on the path (+) for n' = 
1, 2, 3 in Fig. 1. vj - VI = 1/4", vi - = 0.3541/4". Line 
styles and notation are the same as in Fig. 2. 

In Fig. 3 the quadratic corrections in Eq. (1) for 
the correlator ^(?) are zero up to higher orders. The 
quadratic corrections to the first correlator are positive, 
82(1) > 0, and <j)i{t) is convex in In? around /j*. There 



are close-by A2-singularities on the chosen path with 
fq < fq^l = 1 1 2, causing a decay which is encoun- 
tered after the logarithmic decay, which is seen best for 
curves n' = 1 and n' = 2 in the upper panel of Fig. 3 for 
^1 (f)<0.15. While in Fig. 2 the validity of the asymptotic 
approximation around /* was hmited by a higher plateau 
from above, in Fig. 3 the lower plateaus present a bound- 
ary for the application of Eq. (1) from below. Because 
higher plateaus do not interfere, the asymptotic descrip- 
tion appUes three decades earlier in Fig. 3 than in Fig. 2 
for comparable distances from the A4-singularity. 

The estimate for the sign of Ci given above mainly 
depends on the relative distance of the path chosen in 
control-parameter space to the lines Ci = and £2 = 0. 
As the ordering of the latter two Unes and the path is re- 
versed for the case 82(2) = in Fig. 1, Ci is expected 
to be negative there. The path labeled n' in Fig. 1 is now 
closer to the Une £1 = than to £2 = 0. In Fig. 3 we find 
indeed, that a comparison of solutions and approxima- 
tions (1) indicates a negative value for Ci to account for 
a steeper slope {B — Bi—Ci). Despite those small devia- 
tion, the correlators are described well by the asymptotic 
laws. 

The comparison of Figs. 2 and 3 is summarized as fol- 
lows. For each value of q there exists a fine with van- 
ishing quadratic correction B2{q) for the specified q in 
the approximation of Eq. (1). On this line, the logarith- 
mic decay is displayed best for the correlator specified 
by q. Moving to control parameter values above this line, 
Biiq) < 0, introduces concave decay of the correlator 
<l>q{t) in Int. Going below the fine, B2{q) > 0, yields a 
convex decay in the correlator (j>q{t). For increasing the 
value of Kg, the curves specified by B2{q) = rotate 
clockwise around the A4-singularity. 
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FIGURE 4. Glass-transition diagram for the two-component 
model for V3 = 45. Notation and Une styles are the same as in 
Fig. 1. The circle marks the As-singularity. Heavy and light fuU 
lines denote A2-transition points with their respective higher 
and lower than the value f° at the Aj-singularity. The crossing 
point of the two lines is indicated by a square. 



For control-parameter values V3 > V3, the A4- 
singularity is replaced by an As-singularity which 
terminates a glass-glass-transition fine that extends 
further into the arrested region as V3 increases. The 
lines e, = 0, El = 0, £2(1) = 0, and ^2(2) = for the 
A3-singularities emerge from smooth transformations 
of the ones at the A4-singularity. Figure 4 shows a cut 
through the transition diagram for V3 =45. Different 
from the situation in Fig. 1, the line B2(l) = is now 
well separated from liquid-glass-transition lines while 
the line £2(2) = is located completely in the arrested 
regime. The dynamics on the line £2(1) — for this 
A3 -singularity and the related crossing scenario on the 
same Une is discussed in Ref. [13] in great detail. 



TIME SCALE 

It is obvious from Figs. 2 and 3 that the time scales T for 
a given state determined by matching, e.g., (t) at or 
fi + /i can deviate considerably. The resulting scale is 
also different if t is fixed for different correlators. How- 
ever, asymptotically close to the higher-order singular- 
ity, fg approaches zero and also the time scales deter- 
mined by using different correlators converge towards 
each other as explained earlier [13]. Therefore, only the 
scale fixed by (t) = /j* shall be considered in the fol- 
lowing. 

The time scale T when considering only the leading 
order result in Eq. (1) is given by [28] 

log,oT-l/|ei|i/^ (4a) 

for the A3 -singularity and from a similar argument [29] 
for the A4-singularity as 

logioTocl/|£i|l/4. (4b) 

Thus, the logarithm of the time scale is the relevant quan- 
tity to be discussed further on, and logjg T is expected to 
diverge like a power law when the separation from the 
singularity vanishes. Figure 5 shows that this divergence 
is indeed stronger for the A4-singularity than for the A3- 
singularity. For |£i | < lO^-' the time scale T is described 
qualitatively by the asymptotic laws. The laws in Eq. (4) 
were verified quantitatively by extending the analysis 
to times as large as logigT ~ 100. With one exception, 
the time scales deviate to lower values for t for larger 
separations |£i|. Hence, the slope of the logjQ T-versus- 
logioC^^i ) curve is even larger in that regime than given 
by the asymptotic laws in Eq. (4). In the lower panel for 
|£i I > 10~^, the time scales for the A3-singuIarity are de- 
scribed better by the law for the close-by A4-singularity. 

It is seen for the case of the A4-singularity that the 
proximity of a line of Uquid-glass transitions (scales T 
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FIGURE 5. Time scale T determined by matching (pi (f ) at /j* 
or /j° for the paths 82(1) =0 (crosses) and 52(2) = (squares) 
from Figs. 1 and 4, respectively. The dotted lines are guides 
to the eye. The corresponding asymptotic laws (dashed lines), 
T 0= exp[l/|ei I'/"*] for the A4 and T 0= exp[l/|ei | '/^] for the 
A3 are fitted to the data for the paths where 62(1) =0 around 
T 10'"''. The asymptotic law for the A4-singularity is redrawn 
in the lower panel as chain line. 



indicated by crosses) influences also the validity of the 
asymptotic law for t: Only for extremely large times, 
T follows the law (4b). The time scales are also much 
larger on that path as discussed already in connection 
with Fig. 2. 



CONCLUSION 

It was shown that the asymptotic expansion in Eq. (1) 
provides a way to divide the control-parameter space into 
distinct regions by setting the dominant correction to the 
logarithmic decay laws to zero for different correlators. 
If plotted in the form of Figs.l and 4, the glass-transition 
diagrams of the two-component model are easily seen 
to be topologically equivalent to the ones shown for the 
square-well system in Refs. [14, 17]. For increasing Kq 
in Eq. (1), the lines of vanishing quadratic corrections 



rotate clockwise around the higher-order singularity in 
the two-component model, which mimics the behavior 
in the microscopic model for increasing wave vector. 

Using the asymptotic expansion to describe the solu- 
tions of the equations of motion (2), the approximations 
fit rather accurately already in a regime that should be 
accessible to experiments and molecular-dynamics sim- 
ulations. However, the presence of other glass-transition 
singularities in the vicinity of higher-order singularities 
can have drastic influence on the window in time where 
the logarithmic laws are applicable, cf. Figs. 2 and 3. The 
validity of Eq. (1) for ^^(f) around f* can be bound ei- 
ther from above at earlier times, cf. Fig. 2, or from below 
to later times, Fig. 3, depending on the specific location 
of the chosen path in the transition diagram. 

The time scale T asymptotically follows the Vogel- 
Fulcher-type law of Eq. (4). That the relaxation processes 
close to higher-order glass-transition singularities are ex- 
tremely slow is manifested by the fact, that instead of 
T itself, the logarithm of the time scale is the relevant 
quantity that diverges in some power law. The differ- 
ent behavior of T for A4- and A3-singularities is clearly 
seen in the numerical solution of the MCT equations, cf. 
Fig. 5. However, these asymptotic laws only show up 
for extremely long times and small distances from the 
singularity. In addition, other close-by glass-transition 
singularities can significantly modify the validity of the 
asymptotic law especially for further separation from the 
higher-order singularity under consideration. 
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